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Abstract 

Too often microclimate studies in tlie field of cultural heritage are published without any or scarce information on 
sampling design, sensors (type, number, position) and instrument validation. Lacking of this fundamental 
information does not allow an open discussion in the scientific community. This work aims to be an invitation for 
a different approach. 

Three main parameters (temperature, humidity, luminance) were monitored in a selected part of a complex 
construction by an inexpensive self-assembled system along some horizontal and vertical vectors. All data was 
then processed and analyse by chemometric methods. Some measurements of oxygen, carbon monoxide and 
dioxide and pressure were also performed. 

Correlation of some indoor and outdoor data was shown for temperature and humidity. In case of outdoor 
changes the indoor environment reacted with a certain delay which is position-dependent and more evident for 
humidity data. The two observed rooms (Career and Tullianum) behave differently and the hypogean one is less 
influenced by the outdoor environment. Instrument validation before and after the campaign, allows to consider 
detected variations as significant. 

The fundamental importance of Sampling Design and of instrument validation before and after the monitoring 
campaign was enhanced. The choice of two main and two minor vectors allowed detection of different behaviour 
for the two rooms, also permitting to detect for both rooms a trend towards a spontaneous microclimate 
necessary for a conservation project. In the next campaign we will focus on the choice of the best sampling 
frequency to use more sophisticated statistical methods. 



Background 

Microclimate never can be univocally defined in space; 
should one consider the indoor part of a building or 
restrict it to a single room? 

Obviously, correct spatial definition of a microclimate 
should also consider all surfaces (walls, statues, coins, 
etc.); their climatic conditions are indeed influenced by 
that of the surrounding environment, but depend also 
on the material 
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Very often the definition is referred to a building or at 
least to a partially isolated area. 

For a building a mesoclimate' can be defined, which is 
influenced by the building itself, its structure, geometry, 
shape and materials and by the surroundings. 

For our purpose we define a macroclimate, charac- 
terised by locally measured meteor-climatic values, which 
are not influenced by the building itself [1]. 

In the case of historic buildings, one or more centuries 
old, both the meso- and the microclimate should have 
reached an equilibrium condition, that can include small 
natural variation; however higher change and or fluctua- 
tion of climatic parameters might occur due to restora- 
tions, excavations, reinforcements (sometimes with 
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improper materials) or to the presence of humans (visi- 
tors, restorers) making measurement campaigns more 
complex. For conservation such fluctuations mean that, 
for example, a wall structure of stone, ceramics, glass, 
frescos, metals and also organic material which were con- 
served for centuries or millennia due to equilibrium con- 
ditions are put in danger in a very short time range 
compared to the sites life-time. 

The present study concerns with a hypogean, which is 
a very particular part of a building and its investigation 
can be difficult; one reason is that - while identification 
of the spontaneous condition of the research environ- 
ment is particularly important and fundamental for a 
conservation project - it is impossible to define a priori 
for a hypogean the time necessary to reach constancy or 
at least spontaneous periodicity of fluctuation, referring 
to Banham's definitions [2]. 

Once macro- and microclimate are sufficiently defined, 
one can compare both in order to identify their correla- 
tion type, e.g. "Conservative Mode", where the walls or 
the burying soil exclude or strongly slow down any 
exchange; "Selective Mode", where some parameters are 
filtered by the barrier while others are not; "Regenerative 
Mode", where an independent internal microclimate is 
created, which instead leads to retroactive mechanisms 
(typical for modern living environments with high waste 
of energy). For a hypogean the external construction 
often functions as filter for the macroclimate and might 
guarantee good microclimatic and air quality conditions. 
Two significant parameters, for conservation, are not stu- 
died in this research; the humidity distribution inside the 
walls [3] and the presence of microorganisms like fungi, 
bacteria, yeast, moss and also anaerobic organisms [4], 
both related to microclimatic conditions and necessary in 
a conservation project. 

But what are good microclimatic conservation condi- 
tions? Ideal or limit values of microclimate parameters 
are defined in some official Italian publications [5] as 
well as in scientific ones [6], but knowledge on these 
values is still poor. Given that none of these values derive 
from an accurate (long-lasting?) experimental study with 
double blind technique as in the case of new pharmaceu- 
ticals and instead they are often the result of intensive 
expert-arguing [personal communication by one compo- 
nent of Normal commission], there is an urgent need for 
many, accurate, reproducible, comparable works, in order 
to avoid that one day microclimate parameters are 
imposed by law [7]. Unfortunately many papers on 
microclimate studies of indoor environments, especially 
in the field of cultural heritage conservation, could be 
discussed critically, because they do not state the type, 
number or position of the sensors used. Often the Design 
of Experiment (DoE), the result and instrument valida- 
tion is not described. Generally few sensors, sometimes 



only one, are applied and time series obtained, which are 
then analysed by means of control cards and/or univari- 
ate graphs, but surely there is a correlation between tem- 
perature and humidity so a multivariate approach is 
necessary [8]. 

In this preliminary study we propose a, we think, suffi- 
ciently rigorous procedure, but for a too short observa- 
tion period, wherelS inexpensive but accurate data 
loggers were used for monitoring three environmental 
parameters (temperature, humidity, lighting) every 10 
minutes along some vectors in a hypogean constituted 
by two rooms one above the other; chemometry was 
used to extract the most significant information. 

Experimental 

The position of the building in the ancient Forum Roma- 
num is shown in fig.A in additional file 1 (all figures in 
additional file 1 have a capital letter as index); it was 
probably built in the V. or VI. century b.C., above a water 
source. 

The whole building is composed by a vertical series of 
rooms of which we monitored only the lower two, the 
TuUianum and the Career while our study does not 
include the chapel "Santissimo Crocifisso" and the 
church "San Giuseppe dei Falegnami"; complete informa- 
tion on the buildings history, structure and materials can 
be found in a paper of D.B.Kerner [9]. Currently the Car- 
eer constitutes a partially buried church accessible by a 
descending stair with about 20 steps. The TuUianum, 
completely buried, had the floor continuously covered 
with a stratum of water the origin of which goes beyond 




Figure 1 Sensors position according to tine Sampling Design 
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the aims of this study [10]. In fig.l the position of the 
sensors. 

Basing on the archaeologists maps the dimensions of the 
two rooms were obtained; for the TuUianum wall thick- 
ness of about 1.7 m, an internal surface of ca. 25 m^ and a 
volume of ca. 56 m^ were found. For the Career walls 
thickness varies from 1.7 m (external) to 1 m (versus other 
indoor environments), surface is approximately 60 m^ and 
volume ca. 300 m^. 

Luckily, the entire complex is not equipped with a 
heating, ventilation and air-conditioning system (HVAC), 
so, to ensure against anthropogenic variables, it was suffi- 
cient to close both rooms to public for the entire mea- 
surement campaign. The pump used to drain water from 
the TuUianum was set on automatic mode in order to 
avoid that the water level raised up some cm so ensuring 
against sensor diving. 

For microclimate monitoring we used 13 Hobo U12- 
012 by Onset USA for measurements of humidity (RH 
%), temperature (°C) and Luminance (foot-candles or 
lumens/ft^). 

A EB20-THP, by Ebro GmbH, Germany was also used 
as equipped with certified sensors for humidity (RH%), 
temperature (°C) and pressure (mbar). 

In the preliminary coherence study (see paragraph 
Sensors calibration) two VWR cod. 61027-056, certified 
Hg thermometers AA grade, were used for controlling 
and calibrating of data loggers. 

A LabQuest data station with sensors for CO2, O2 and 
temperature measurements was used as well (C02-bta, 
02-bta and STS-bta by Vernier Tech., USA). 

Two Lascar data loggers (UK) were used for testing 
(EL-usb-CO for CO in air and EL-usb2 for RH%, °C). 

In order to compare micro- and macroclimate, 
meteorological (macroclimate) data was obtained from 
"Osservatorio Meteorologico del CoUegio Romano", 
active since 1782 and located in the centre of ancient 
Rome and on the top of a 63m (above sea level) tower, 
situated in 41.900°N, 12.480°E. TuUianum instead is situ- 
ated in 41.893°N, 12.484°E. 

Data elaboration software were mainly free ones, e.g. 
Past, Datalab, Gnumeric, WinlDAMS, Libre Office and 
some commercial ones like MVSP and Lotus 123. 

To read data from data-logger producer software were 
used (available on brand web site). Only one software 
was not free and we apologize for this. 



Methods 

Our work focused on the monitoring of two principle 
parameters, temperature ("C) and relative humidity (RH%). 

As above said, an other parameter strongly influencing 
the conservation state of an object is bound to lighting. 
In a running research we are considering the effect of 
both the spectra of the light source and its wavelength- 
weighted power while in the present study only the last 
was monitored; illuminance was measured due to its 
greatest significance in studies regarding frescoes, paint- 
ings and any object in general. 

15 data loggers were employed as passive samplers to 
monitor the three above said parameters; they were 
easily fixed also along vertical vectors (weight 46g) with- 
out disturbing the spontaneous environmental air flux. 

When passive samplers are employed sampling fre- 
quency must be lower than the sensors rising time, so 
we choose a 10 min sampling time (see next paragraph) 
as the temperature rising time for Hobo is 6 minutes, 
see table 1. 

Even if active samplers may be used at a higher fre- 
quency with respect to passive ones, we choose the last 
as, in our opinion, for a long monitoring period the 
adopted frequency allows to obtain significant informa- 
tion and, overall, the artificial air flux of known flow 
rate created by a pump and oriented on the active sen- 
sor [11], involves a self-heating especially in restricted 
indoor environments. Really, the self heating of a ther- 
mistor is a problem well known to electrical engineers, 
even a small 2W pump motor can modify the microcli- 
mate and thus measured values, when employed for 
longer periods. 

In this first campaign an experimental data station 
with some sensors was also used in order to measure 
gas concentrations (e.g. CO2 and O2) in air. Further an 
attempt was made to evaluate the CO level as an indica- 
tor for contamination of the hypogean with external, 
urban, air. In order to identify the "spontaneous" values, 
the two environments were not entered during the 
15 days measurement campaign. 

Sampling design 

In this work the choice of vectors direction and position 
for sensor placement was based on Sampling Design 
(SD) using an expert-sampling among all authors for 
vector position; on each vector a systematic-sampling 



Table 1 Comparison of the values of the 3 data loggers, according to the producer 

hobo lascar ebro 

parameter temp RH% light intens. temp RH% temp RH% mbar 

accuracy +/- 0.35 +/- 2.5 - +/- 1 +/- 3.5 +/- 0.5 +/- 3 +/- 5 

resolution 0.03 0.03 0.1 0.5 0.5 0.1 0.1 1 

drift, one years 0.1 1 - n.a. 1 0.5 1 5 

rise time, min 6 1 - 0.2 0.1 5 2 1 
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with even spaced distances including the extremes, was 
applied. 

Position of the two main vectors is vertical (trespas- 
sing the two rooms, 6 sensor points) and horizontal (in 
the basement, 6 sensor points); the minor two vectors 
are one along a horizontal axis in the upper room (3 
sensor points) from the door to the altar and one in the 
basement (3 sensor points) from the stairs to the oppo- 
site wall (fig.l). 

The long horizontal vector is designed at about 50 
mm from the floor and at the cross of the two vectors, 
the concentrations of CO, CO2 and O2 were measured 
too. 

In the time domain we used a systematic-sampling 
with a reading every 10 minutes, but the choice of the 
proper sampling frequency was based on a preliminary 
24h measurement campaign, using the faster available 
sensor, at 2 readings each minute. 

Obviously immediately after the campaign one notices 
that some additional data logger somewhere would have 
been useful, leaving space for further campaigns and 
making work even more complex. 

Sensors calibration 

Sensors control before and after measurement cam- 
paigns is a minimum requirement which is generally not 
stated in the microclimate studies on Cultural Heritage 
Sites. 

Especially when inexpensive sensors are used problems 
could arise due to lacking of coherence of the single mea- 
surements; the same data-loggers are often available in a 
certified version, but as our aim was also the assembling 
of a inexpensive measurement system, the price of these 
latter was obviously too high. So, at the same cost we pre- 
ferred a double number of data points buying more data- 
loggers and to verify by our own the sensors quality. 

To do this, in a preliminary measurement campaign 
all 15 sensors where placed for a week in a hypogean 
ambient similar to the one to be studied, near to each 
other and around two certified instruments (the EB20- 
THP and VWR Hg thermometer); in this way we 
checked the sensors rising time, accuracy, hysteresis, 
and coherence before starting our measurement cam- 
paign. This allows to estimate the sensors response to 
environmental variations, e.g. a coherence of 0.5 °C 
allows us to regard a variation of 0.8 °C between two 
rooms or along a vertical axis as significant. 

Once the sensors coherence was confirmed they were 
reprogrammed and positioned in the Tullianum for our 
study purpose. 

At the end of our work all sensors were re-checked in 
the same way and in the same place and results are 
shown in additional file 1 as fig.B (temperature) and fig.C 
(humidity); the graphs show that no damages or 



calibration lacks occurred to the sensors during the mea- 
surement campaign because data are coherent consider- 
ing the accuracy limit of the data-loggers (table 1). A 
Vernier data station was also used equipped with tem- 
perature, CO2 and O2 sensors (fig.D in additional file 1, 
CO trend is also reported); temperature data showed an 
offset of ca.0.5 °C and was therefore excluded even if 
used by the data station to apply the right correction fac- 
tor to data coming from the two other sensors. 

Data treatment 

The obtained data was put in a 3D-matrix composed by 
16 columns, with sensors as variables, 2150 rows for 
each sampling time and 7 layers for temperature, RH, 
luminance, barometric pressure, CO2, O2, CO. Some 
layer columns are empty as not all data loggers can 
detect all parameters. For example the barometric pres- 
sure layer has only one column (Ebro sensor). 

The above described matrix was then studied by 
Exploratory Data Analysis in order to identify any corre- 
lation of indoor and outdoor measured values as well as 
to identify the rooms spontaneous values. 

Data of each measured parameter (layer) was plotted 
in monovariate graphs. Normally Run-Sequence plots of 
parameters (temperature, humidity, light intensity..) ver- 
sus time or a function of time (here minutes from mid- 
night of the first day). 

In these graphs prefixed limits can be added and She- 
whart Charts are obtained (also known as Control 
Charts), which enhance the environments characteristics. 

The successive step in EDA is data analysis by Box & 
Whisker Notched Box Plots of a single parameter (a 
matrix layer), where variables (sensors) are plotted on 
the X-axis and raw measured distribution on the Y-axis. 

This powerful analytical tool shows, here in the 
McGill [12] version, in the median domain, the distribu- 
tion, spread, skew and percentile values 5% and 95% by 
whiskers; no outliers are shown. In our case only one 
matrix layer was plotted at a time so avoiding problems 
which can arise when variables with different scaling are 
plotted together. 

Successively the three Scatter Plot Matrixes were 
designed showing up all possible Var-Var graphs 
between all variables (sensors), starting from raw data, 
in order to enhance correlations as well as distribution 
type of each variable on the diagonal. 

On the diagonal of the Scatter Plot Matrix a histo- 
gram of ten bars is shown for each variable, which 
roughly allows identification of the distribution form; 
the values on histogram show the height of the major 
bar as the total is normalised to 1. 

This important graph is often substituted by a correla- 
tion matrix, using the r^ Pearson coefficient; this latter 
works only on linear data, which is not always the case. 
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especially in the cultural heritage sector (please see the 
complex column of HobolO, on the altar, on the next 
figure). 

In temperature and humidity layers there is also a col- 
umn with macroclimate (outdoor) data; correlation of data 
in this column and those of the others is possible and 
plausible. However macroclimate data has a 30 minutes 
sampling interval; therefore a new matrix was constructed 
in which the indoor data was downsampled in order to 
allow better comparison of indoor and outdoor data. 

Correlation of variable couples was then analysed and 
this was done not only by computing the common r^ 
coefficient, which may not be sufficient. An example of 
calculated regression parameters can be seen in table 2 
for barometric pressure. 

If one plots in graphs parameters which best describe 
correlation, in terms of quality, here rows enhanced in 
blue, it is possible to compare two or more correlation, 
which differ for instance for the data amount (many 
validation techniques like LeaveOneOut or Bootstrap 
are based on these descriptors). 

A building isolation can be evaluated also by Loss/ 
Gain for Transmission graphs, where internal and exter- 
nal data are shown and compared to a perfect correla- 
tion which is shown as graph diagonal. Any sign above 
the diagonal is a gain of the building compared to the 
outdoor environment, e.g. can be due to a heating sys- 
tem for a temperature graph. By using labels for single 
data points these graphs allow also evaluation of hyster- 
esis, spontaneous behaviour, changes with time. The 
sensors accuracy is fundamental in these graphs. 

A temperature-humidity correlation study was the next 
analysis step. Thus two columns of two overlaying layers 
were compared and for each point variation of both para- 
meters as well as the 98 % confidence ellipses were 
obtained. 

All these data exploration methods together allow eva- 
luation of the buildings isolation from the outdoor 
environment, as well as enhancing internal gradients; 



further spontaneous values, not normal behaviour and 
sudden changes, which are the most dangerous events 
for a cultural heritage site, can be evaluated. 

Results 

Time plot 

Microclimate monitoring results are represented, as run- 
sequence-plots, in fig.2 (temperature in degree Celsius) 
and fig.3 (relative humidity percentage). In order to 
allow better data overlapping, macroclimate data (red 
dotted line) were added in the same graphs (red scale 
on the right). 
These simple graphs already enhance: 

- a more or less evident correlation between indoor 
and outdoor (macroclimate) data. 

- a certain delay, more evident for humidity data, for 
indoor changes with respect to outdoor changes. 

- as expected, the two rooms behave differently; the 
hypogean is less influenced by macroclimate and show a 
lower variability for both the measured parameters. 

- a tendency versus a spontaneous microclimate, more 
evident for humidity so putting in evidence the need to 
extend in time the next measurement campaign. 

- finally, the coherence among sensors was checked 
before and after the campaign, (examples of the latter in 
fig.B and fig.C in additional file 1), allowing to consider 
all variations shown in graphs of fig.2 and fig.3 as 
significant. 

- careful graph observation allows stating that for 
these two different rooms the attempt to draw univari- 
ate Shewhart or Control Charts and define limit values 
to be respected by a HVAC would be confusing and 
stupid. 

In order to improve data reading one can draw a ser- 
ies of graphs each one relative to a single day. As an 
example, in fig.4(a,b) temperature and humidity data, 
along the vertical vector for day 4 and 5 respectively, 
enhance the parameters different behaviour according to 
Benham [2] and the above said delay. 



Table 2 Evaluation of the delay between Internal and external barometric pressure studying the lag correlation. 



parameters \ lag (min) 


0 


30 


60 


90 


120 


150 


180 


sample covariance 


13.644 


13.631 


13.593 


13.547 


13.496 


13.437 


13.369 


population covariance 


13.625 


13.612 


13.574 


13.528 


13.477 


13.419 


13.350 


Pearson product-moment correlation coefficient, sample 


0.990 


0.990 


0.988 


0.986 


0.983 


0.980 


0.976 


square of the Pearson correlation coefficient 


0.981 


0.980 


0.977 


0.972 


0.966 


0.960 


0.952 


assess of Std Err of Pearson's r 


0.005 


0.005 


0.006 


0.006 


0.007 


0.008 


0.008 


variance of post-fit residue 


0.259 


0.270 


0.317 


0.380 


0.451 


0.536 


0.640 


residue sum of squares (RSS), sample 


183.467 


191.739 


224.471 


269.190 


319.809 


380.278 


453.889 


root mean square deviation (RMSD), sample 


0.508 


0.520 


0.562 


0.616 


0.671 


0.732 


0.800 


standard error of the Y estimate 


0.509 


0.520 


0.563 


0.616 


0.672 


0.732 


0.800 


computed Slope 


0.972 


0.971 


0.968 


0.965 


0.961 


0.957 


0.952 


computed Constant 


27.831 


28.801 


31.486 


34.821 


38.492 


42.711 


47.662 
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Figure 2 Trend of the temperature during the measurement campaign by all the sensors and comparison with macroclima 



Study of the values distribution 

For remarking of each distribution and mutual beha- 
viour Box& Whisker plots can be used starting from raw 
data. 

Including all sampling points the graph in fig.Sa (tem- 
perature) is obtained where two main groups can be 
observed, each one including the sensors located in each 
of the two different room. Really, sensors 1, 2, 4 and 10, 
located in the upper room, behave completely different; 
wider ranging values and higher absolute values for the 
medians in comparison with those located in the TuUia- 
num (ebro, lascar, 3,5,7,8,9,11,12,13) are clearly evi- 
denced. Sensor 6 is also located in the Career but just 
on the grid connecting the two rooms and it is not so 
easy to insert it in one of the two groups. More in this 
group the sensors which are part of the horizontal vec- 
tor (1, 4, 10) behave similarly for what concerns both 
the median value and range while on the vertical vector 
the expected trend can be observed; really, coming 
down, the graph shows a decreasing trend for both 
median values and spread as well as decreasing differ- 
ences between points. Along the main horizontal vector 
of TuUianum a very similar behaviour is evidenced with 
the exception of point 7 that shows a slightly lower 



spread and higher median while on the minor vector 
temperature is almost constant but ranges slightly 
different. 

In fig.Sb, the humidity distribution evidences the same 
two main data groups well separated; one characterised by 
wider ranging values and a lower median (points 1, 2, 4 
and 10) and another one including all remaining data 
points. Sensor 6, in this case, is clearly included in the Tul- 
lianum group even if, together with sensor 5, located on 
the stairs connecting the rooms at about 2 m from the 
lower floor, shows a similar median value enough low 
with respect to all the others. The presence of the sub- 
group including sensors 11, 12 and 13 is also better evi- 
denced; such sensors are positioned on the main 
horizontal vector, inside a tunnel probably coming versus 
the Cloaca Maxima and connected to the room through a 
metallic door that was left open during the campaign. 

Correlation among points 

As described in methods section. Scatter plot matrix were 
designed in order to detect possible correlations among 
the columns of the single matrix layers. Starting with lumi- 
nance data in fig.6 one observes that only sensors n° 1, 2, 4 
and 10 show a series of values (see columns I, II, IV and 
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Figure 3 Trend of the relative humidity during the measurement campaign by all the sensors and comparison with macroclima 



X); values of sensor n° 4 are distributed in a very narrow 
range and point 6 has only two values, confirming that the 
used data loggers do not have a LOD sufficient for the 
luminance measurements due to the scarce photon num- 
ber in this hypogean. 



Fig, 7 shows the matrix of temperature plots and it 
enhances immediately the low resolution of the Lascar 
sensor (very similar values compared to those measured 
by the other sensors). All other sensors show a more or 
less evident positive correlation. Sensor n° 8 for example, 
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Figure 5 Box&Whisker plots for temperature (a) and relative humidity (b) measured by all the sensors 



situated below a grating connecting the Career and the 
Tullanium, shows a high hnear correlation (r^ = 0.9819) 
with point n° 6 and all TuUianum points; a significatively 
lower correlation with Career points is evidenced (r^ = 
0704 with point n° 2). 



In fig.8 the scatter plot matrix for relative humidity 
is shown; even if less evident with respect to fig.7, a 
positive correlation among sensors can be observed as 
well as a wider variability for the entire matrix. It can 
be noted that, as expected, sensors 11, 12 and 13, 
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Figure 6 Scatter plot Matrix for luminance measured by all the sensors 
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Figure 7 Scatter plot Matrix for temperature measured by all the sensors 



located in the tunnel and on the Tullianum main vec- 
tor are well correlated with each other and also with 
sensors 9 and Lascar that are located on the same vec- 
tor; even if sensor 3 is also located on the same vector 



and closer to the tunnel, a lower correlation is evident 
that can be attributed to a position close to the stair 
that allow a higher air flow as already evidenced 
above. 
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Analysis of the structure's isolation from the external 
environment cannot be performed by a simple scatter- 
plot, instead one has to calculate all correlation para- 
meters [13] among two columns of the same layer. 

Daily variation 

Before starting lag- correlation analysis it seems useful to 
identify time-windows hidden in time-series (repetition 
of the same values with a certain frequency due to cyclic 
phenomena such as night-day or seasonal variations). A 
number of software allow detection, expensive ones and 
free ones like Past [14]. 

Time series analysis are performed in the same way as 
Spectral Analysis, starting with Fourier Transformation 
which allows shifting from the time domain to the fre- 
quency domain. For temperature and RH data series of 
sensors 2 (top of Career) and 8 (top of Tullianum), a 
time window of 1430 minutes is detected for the first 
peak of the FFT spectra. 

Usually in combination with spectral analysis Sinusoi- 
dal Regression (often called Sinusoidal Curve Fitting, 
SCF) is used, where a sum of n sinusoids with specified 
periods is used to fit a known list of time dependent 
data. Past algorithm of SCF is based on a least-squares 
criterion and singular value decomposition. 

This can be useful to check and model periodic phe- 
nomena in time series, as in our case. Further it is inter- 
esting to compare the results of FFT and SCF in order 
to obtain result validation and to avoid overfitting. 

Expressing the identified time window in hours, one 
finds a values close to 24 (1440 minutes), which con- 
firms influence of the external environment on our 
indoor ones. 

To check fast change in values (temp. RH, lux) we 
must use a sampling frequency more than two times 
faster than the change under study (Nyquist-Shannon 
sampling theorem), but this is difficult with a sensor 
with rise time of 6 or 10 minutes, one other big pro- 
blem in microclimate monitoring often silenced. 

Lag correlation 

The meteorological data used had a 30 minutes sampling 
interval, thus our measured microclimate data had to be 
reduced by downsampling as described in methods 
section. 

For correlation analysis of meteorological data (out- 
door, X axis) and measured microclimate data (indoor, 
Y axis) Ordinary Least Square was used; data of the Y 
axis was then shifted by 0, 30, 60, 90 and 120 minutes 
and obtained parameters registered in new matrixes, like 
the one in table 2. 

Correlation coefficients were calculated for all avail- 
able data without using a time window in order to avoid 



any anomaly. If number of data is much higher (10-30 
times) than the width of time window this techniques is 
sufficiently robust to avoid dependence from some not 
normal data (e.g. presence of a small heater on day 12 
or a window open all night on day 18). 

In contrast to other works our first approach was data 
validation and results coherence; correlation analysis of 
barometric pressure (mbar) of Tullianum data and 
meteorological data was performed first. As no HVAC 
is present in the Tullianum, barometric pressures must 
be equal, see fig.F in additional file 1. 

In table 2 obtained coefficients are shown. A value of 
r^=0.9665 is observed, which is very good considering 
the distance of the two measurement points (850m) and 
the different instrument types and cost. 

Further, lag-correlation is maximum for t=0 while all 
blue values worsen the situation, introducing a delay on 
the values of the Y axis, which are those measured in 
the Tullianum. 

Computing all correlation parameters listed in table 2 
for all space points, for temperature and humidity data, 
graphs containing the data referring to all vectors can be 
draft. 

As an example figs.9(a-d) show a series of graphs 
obtained from Table 3 (in additional file 2), using para- 
meters marked in blue, for temperature and humidity 
measured by sensor 8, situated below the communicat- 
ing hole of the two rooms. All parameters evidence a 
delay for temperature equal to 90 minutes, but with a 
determination coefficient r^=0.224; thus only 22.4 % of 
the variance of Y can be explained by X, this means that 
the system must be considered an almost completely 
isolated one. 

Humidity always shows the same delay with a 
r^=0.061, thus X and Y must be considered independent. 

A simple visual comparison is not sufficient and cer- 
tainly the common analysis of delays between minima 
and/or maxima (if detectable) is not useful in order to 
obtain a correlogramm. 

The lag correlation is already evident in figs.4a and 4b, 
i.e. the line plots of the raw data obtained with a 1440 
minutes time window for temperature data on day 5 and 
for RH data on day 4 respectively (days 4 and 5 were cho- 
sen due to strongly varying values of parameters). 

The next step in our analysis on the building was to 
show up r^ values and slopes of the linear correlation of 
each main axe obtained by comparing outdoor and our 
indoor sampling points using a proper delay prior iden- 
tified as above described. 

In fig. 10a, trends for r^ and slope values, relative to 
the correlation between macro- and micro-climate on 
the vertical vector are shown. For both temperature and 
humidity, it is well evident a decrease starting from 
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Figure 9 Trend of the correlation parameters as a function of the delay between measures performed by sensor 8 and macroclima: a) 

Pearson's coefficient of determination, r^; b) variance of post-fit residue; c) residues sum of squares (RSS); d) root mean square deviations on 
residue (RMSD). 



position of sensor 2 for both parameters and, as 
expected, starting values significantly higher as well as a 
quicker decrease for sensors located in the Career. 

Fig. 10b shows lower values for both r^ and slope and a 
less evident trend for both parameters relative to the 
series of the shorter horizontal Career vector. 

A flat trend and a further decrease for both slope 
(close to zero) and r^ values can be seen in fig. 10c and 
d relative to the two horizontal TuUianum vectors, 
remarking that the ambient has to be considered an iso- 
lated one. 

For all the vectors, lower slope and r^ values are evi- 
denced for humidity with respect to temperature; such 
behaviour is consistent with a greater influence of the 
Tullianum humidity with respect to the macroclima. 

Insulation 

Rooms isolation can be estimated as well by a Loss/Gain 
representation; firstly a daily mean is calculated for both 
the outdoor (meteorological) and the indoor data (the 
latter for both the Career and the Tullianum separately, 
considering all the relative sensors) then correlation 
graphs were drawn for temperature and humidity data. 



Fig. 11a and lib shows the temperature and humidity 
graphs. On temperature, with few exceptions (as an 
example Career on days 19, 29 and 23), both the Tullia- 
num and the Career show lower values compared to the 
external ones. Similar daily variations, enhancing com- 
munication between the two rooms, can be seen for 
Career and Tullianum; more, lower variations can be 
seen in the order: outdoor> Career > Tullianum (i.e. maxi- 
mum variation result about 5, 2 and 1 °C respectively), 
enhancing the expected higher isolation for the 
Tullianum. 

To understand at the best the above behaviour, it is 
useful to consider that the position of the building 
entrance, that is the only point in direct communication 
with outdoor through an open-work door, is east- south- 
east oriented, so, sunlight can enter from about 8 until 
about 14. Further, it must be pointed out that Tullia- 
num is a classical hypogean but also the Career is par- 
tially underground; as a consequence, illumination level 
is low everywhere as can be seen in fig.E in additional 
file 1 where only the upper sensor on the vertical vector 
(Hobo2) seems to be stimulated. On the other hand the 
underground position of the Career allows a lower 
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Figure 10 Slope and trend for the best correlation among delay relative to: a) vertical vector; b) Career horizontal vector; c) Tullianum 
main horizontal vector; d) Tullianum minor horizontal vector. 



influence of eventual external wind through the open- 
work door. Finally, all the walls are very thick so ensur- 
ing a good insulation. 

An expected opposite trend, i.e. value are above the 
outdoor ones, is shown in fig. lib for humidity with the 
exception of data relative to Career on days 17, 21, 23, 



24 and 29. The same flatting trend on variations, above 
seen for temperature, is evidenced with maximum value 
of about 25, 18 and 10 % for outdoor. Career and Tul- 
lianum respectively. Really the presence of groundwater 
ascending in the Tullianum and the huge walls funded 
directly on the water layer cause high internal humidity 
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levels, like the one shown in the final part of graph in 
fig.3; any reasonable conservation concept should not 
interfere on this. 

Multivariate 

The reply at the question "is there a correlation of 
humidity and temperature?" is obviously "yes" but it is 
not unique and exhaustive. A first answer can be based 
on fig. 12, showing a scatterplot, after column centering 
of the outdoor data (humidity versus temperature) dur- 
ing our measurement campaign; a not so clear reply can 
be given considering the indoor data. 

Column centering is an important data pretreatment 
because it leads to scatterplots which show directly the 
swing of the values on the axes; for fig. 12 we obtain 
about 14 °C for temperature and about 70% for humid- 
ity with a negative correlation. Further outliers, trends 
and leverage are easily enhanced while centering does 
not modify distribution structure; instead it can lead to 
different results in the successive EDA phase [15]. 

In order to enhance the study of the distribution form 
for temperature and humidity data for each sensor, one 
can compute bar charts with class distribution. Choice of 
class width is not a trivial task, should consider values well 
above the instruments resolution level and producing a 
number of classes ranging from 10 to 30. Here we chose a 
width of 0.5 °C and of 2 % for temperature and humidity 
respectively, but different choices can change "the shape". 

Multivariate representation like that in fig. 13 for Career 
and fig. 14 for Tullianum was calculated for each measure- 
ment point and also the class distribution was draft. In 
figs. 15a and 15b the distribution of graphs for the sensor n 
° 1 (behind the door, see fig.l) are shown. Confidence 











: 



- \ t u..^ . ^i. \ 











-12 -9 -6 -3 0 3 6 9 12 



t-meteo 

Figure 12 Correlation graph between outdoor temperature and 
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ellipsis of fig. 13a (and of all the other ones here reported) 
was calculated for 98%; it shows an enough wide variation 
of both parameters (spread equal to 5.4 °C and 24 %RH) 
and no clear correlation that anyway seems to be more 
positive than negative; class distribution of temperature, 
shown in fig. 15a results enough regular with a maximum 
around 15 °C while in the humidity graph (fig. 15b) higher 
value are significatively prevalent. 

The complete series of figures are available in addi- 
tional file 1, from fig.G-a to fig.V-c. 

Among this huge amount of data, one should focus on 
data from sensor 2 (fig. 13b), located at the top of the 
vertical vector (see fig.l); a wider ranging value for both 
parameters is evident with respect to the previous sen- 
sor (spread equal to 6.4 °C and 37%RH) and a no clear 
correlation could be seen. 

Interesting is also point 5 (located in the lower part of 
the inner stair, see fig. 14a; this figure shows a narrower 
and high sloping ellipsis with respect to the previous two; 
it must be also noted an enough crowded queue, out of 
the ellipsis, where a humidity variation of about 5% corre- 
sponds to an almost constant temperature. Humidity 
spread results equal to 20% so being no so different from 
that of sensor 1; an other similarity with such sensor is the 
positive correlation while temperature spread results sig- 
nificatively lower and equal to 2.3 °C. Looking for a simi- 
larity due to the position, it can be noted that, with respect 
to Career, sensor 1 is the nearest to outdoor and sensor 
5 is the nearest to Tullianum; it is evident that macroclima 
have the maximum influence on the first while the effect 
of the very humid Tullianum is higher on the last. In 
fig. 13c, relative to sensor 6, one can observe a near spheri- 
cal shape for the ellipsis and a queue, similar to the one 
discussed above; in fig. 14b, relative to sensor 8, a narrower 
shape along y axis results from the grouping of points 
inside the ellipsis and the external queue is more similar 
to the one in fig.l3c. The two sensors are located above 
and below the hole connecting the two rooms and the 
above trends evidence the influence of the high humidity 
inside Tullianum on that of the Career. Basing on fig. 14c 
it can be stated that, as expected, no influence from the 
macroclima is evidenced by sensor 11, located in the far est 
position inside the tunnel toward Cloaca Maxima; in parti- 
cular, humidity is practically 100% for most of the points 
(see fig.R-c in additional file 1). 

The class distributions for all the sensors (fig.G to fig. 
V, letters b and c in additional file 1) show a decreasing 
number of classes in the order meteo> Career > Tullianum 
with decreasing temperature and increasing humidity in 
the same order. 

As well known and already said in the introduction, 
major suffering of any ambient is due to the rapid 
changes of microclimatic parameters which cause a ser- 
ies of chemical and physical processes in the structures 
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Figure 13 Correlation graphi between indoor temperature and relative humidity measured by: a) Hobo 1; b) Hobo 2; c) Hobo 6; d) Hobo 
5; e) Hobo 8; f) Hobo 1 1 



as well as in mobile objects; instead, all slow changes are 
certainly sustained by materials which survived in these 
environments for 300 or 2000 years. In order to detect 
eventual rapid changes FFT and SCF can be very useful. 
An example is given in fig. 16 (humidity from sensor 6 
with 5 min. sampling interval) where some enough 
rapid changes are evidenced. Such changes can be seen 
better using the FFT, in fig. 17 the already identified 
daily change is well evident but also some less frequent 
variations at both lower and higher frequency result sig- 
nificative (red line refers to 95% confidence). 

Conclusions 

Referring to the Tullianum the main exchange with out- 
door seems to take place along the stairs and not 
through the hole in the ceiling. Really figs.5(a-b) evi- 
dence significative differences between sensors 6 and 8 
that are located, at only about half meter each, above 
and behind the grid connecting the two rooms; this can 
be explained with a scarce air flow through the grid 
and/or with a bigger influence of the Tullianum micro- 
clima on sensor 6. An air flow from the tunnel coming 



toward Cloaca Maxima, together with a water layer 
often present on the floor, seem to be the main respon- 
sible for the Tullianum microclima; anyway, also an air 
flow from the stairs plays a role as humidity never 
reaches 100% that, on the contrary, is the almost con- 
stant value evidenced in the tunnel. More, the low 
macroclima influence is also evident in figs.2 and 3. 

Taking into account that the building was always 
closed during the measurement campaign, figs.2 and 3 
evidence a trend coming toward spontaneous conditions 
with lower and at lower frequency variations for both 
temperature and humidity, in fig. 18 humidity data regis- 
tered from midnight of 26 April until the end of the 
measurement campaign evidence at the best what said 
above. Matrix layer was subjected to row-centering to 
enhance the difference among sensors. 

In fig.l8 the same trend is evidenced for the Career 
except for sensor 2 that, as already seen, is located in 
the position most influenced by the macroclima. Such 
figure also reveals that spontaneous values flat toward 
values more closer to those of Tullianum probably as a 
results of the good insulation. 
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Figure 14 Distribution of microclima data measured by Hobo 1: a) temperature; b) relative humidity 
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But the mam question of a conservator is: such 
conditions are safe for the artwork? 

We cannot give the right reply! Really this paper can give 
an input to perform a monitoring at the best but time 
was too short and, as a consequence, no seasonal but 
only daily variations were obtained; further, it was per- 
formed not in real conditions that foresee the entrance of 
visitors with all the consequent problems. Probably a 
longer monitoring would allow to obtain more accurate 
values for the spontaneous conditions an surely in such 
conditions the bigger damage due to variations can be 
eliminated, but who can ensure that spontaneous condi- 
tions would be the best for all the materials present in 
the building? Really in the building are present different 
tuffs, frescoes, stone column each material has 



different chemical-physical characteristics so being differ- 
ently influenced by the environment. Lastly but not in 
importance: what conditions can ensure against damage 
from biological source? The presence of visitors allows 
microorganism entering the building, increase the carbon 
dioxide concentration and so on, but the absence of visi- 
tors cannot ensure that no biological damage can occur. 
More, in our case we have a sure capillary rising that can 
allow presence of salts that must be taken into account. 

Probably is exact the statement of F.G. Ojea is exact "la 
conservacion de las pinturas exige un mantenimiento de 
la temperatura y humedad a niveles lo mas constantes 
posiblesy y de un orden aproximado al que tenia el monu- 
mento antes de ser abierto, condiciones optimas para su 
conservacion, y que, obviamente, sifueran desfavorables 
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no permitinan que el monumento llegara hasta nosotros" 
[16] and it is applicable to our hypogeum. 

What must be done? A correct conservation project 
foresees a perfect knowledge of all the materials that 
must be saved and of all the surrounding to find a com- 
promise between the best environmental conditions for 
all. In our case, surely frescoes are the most precious 
artworks and also those more sensitive to damage; 



surely the microenvironment must take into account 
this. A compromise must be adopted to allow visitors to 
enter the building but surely microclima cannot be pro- 
jected for their well-being, as stated in some book. 

In practice the right procedure is running as the Sov- 
rintendence asked for cooperation to "La Sapienza" Uni- 
versity to carry out analysis of constituent materials, 
map of salts and microclimate monitoring. 
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For the last a new campaign is running in the same 
seasonal period to evaluate the influence of visitors. 
Other rooms will be monitored to evaluate the recipro- 
cal influence and seasonal campaign will follow. 

Also citing by reference [17] '\.,not only is the rela- 
tionship among all the variables at any one time impor- 
tant, but so is the entire past history of the trajectories of 
all these variables,'' confirm the use of multivariate ana- 
lysis as necessary [8]. 

We may conclude this work stating that careful con- 
servation of an entire ambient, here hypogean, declared 
cultural heritage, should foresee a triple action: 

a) identification of real conditions microclimate over a 
long time-scale, with season changes and lasting possibly 
years; 

b) comparison of the identified conditions with those 
indicated in official publications (see [18] for a list) in 
order to plan possible site tailored interventions, prefer- 
ring passive methods; 

c) publication of the done work in international jour- 
nals in order to allow wide knowledge diffusion; work 
description should include details on measurement pro- 
cedure AND RAW DATA to allow eventual further 
data elaboration by of other researchers. 

Additional material 



Acknowledgements 

This article has been published as part of Chemistry Central Journal Volume 6 
Supplement 2, 2012: Proceedings of CMA4CH 2010: Application of 
Multivariate Analysis and Chemometry to Cultural Heritage and Environment. 
The full contents of the supplement are available online at http://journal. 
chemistrycentral.eom/supplements/6/S2. 

Author details 

^Rome University, La Sapienza, Math. Phy. Nat. Science Faculty, Piazzale A. 
Moro 5, 00185 Rome, Italy. '^Cultural Heritage Superintendence of Rome City 
Council, Pz. Lovatelli 35, 00186 Rome, Italy. 

Authors' contributions 

Among the authors, who contributed in the same manner to this work, we 
would like to differentiate S.D.G., our graduating student, whose thesis is 
focused on the chemical/physical characterisation of some Tullianum parts; 
further P.F. is responsible for the cultural heritage site for the Roman City 
Council, who contributed with information on the historic/artistic aspect. 
The remaining authors contributed in a different manner but at the same 
level to this work. 



Competing interests 

None of the authors received any financing, scholarship, contract or 
consulting from societies mentioned in this paper. Instruments were bought 
regularly through national research funding. 

Published: 2 May 2012 
References 

1. Nilsson LO: Interaction between microclimate and concrete a 
prerequisite for deterioration. Constr Build Mater 1996, 10(5):301-308. 

2. Banham R: The Architecture of the Well-tempered Environment. 
Architectural Press, London;13: 978-0226036984 1969. 

3. Hedenblad G: Moisture permeability of mature concrete, cement mortar 
and cement paste. PhD dissertation Lund University, Sweden; 1993, http:// 
www.dissertations.se/dissertation/aalef97c3b/accessed at feb 2011. 

4. Cataldo R, Leucci G, Siviero S, Pagiotti R, Angelini P: Diagnostic of the 
conservation state in the crypt of the Abbey of Montecorona: biological, 
microclimatic and geophysical evaluations. J Geophys Eng 2009, 
6(3):205-220. 

5. . , List of the UNI-Normal available at: http://webstore.uni.com/unistore/ 
public/searchproducts?action=search&usecache=false&committee=BENI 
+CULTURALI+-+NORMAL, accessed Feb 2011.. 

6. Delogu S: Qualita microclimatica e consumi energetici nei musei. Thesis 
dissertation Politecnico di Torino, Italy; 2010. 

7. Strada M, Carbonari A, Peron F, Porciani L, Romagnoni P: The microclimate 
analysis of tezone '105' of the Venetian Arsenale. J Cult Merit 2002, 
3(1):89-92. 

8. Garcia-Diego F-J, Zarzo M: Microclimate monitoring by multivariate 
statistical control: The renaissance frescoes of the Cathedral of Valencia 
(Spain). J Cult Merit 2010, 1 1 (3):339-344. 

9. Karner DB, Lombardi L, Marra F, Fortini P, Renne PR: Age of Ancient 
Monuments by Means of Building Stone Provenance: a Case Study of 
the Tullianum, Rome, Italy. J Archaeol Sci 2001, 28(4):387-393. 

10. Fortini P: Career tullianum il carcere mamertino al Foro romano. 
Mondadori Electa;8843566237 1998. 

11. Ryhl-Svendsen M: Luftschadestoffe in Museen, Eine Einfuhrung in 
Wirkungsweise, Monitoring und Kontrolle. Restauro 2001, 8:613-619. 

12. McGill R, Tukey JW, Larsen WA: Variations of Box Plots. Am Stat 1978, 
32(1):12-16. 

13. Pearson K: On lines and planes of closest fit to systems of points in 
space. Philos Mag 1901, 2(6):559-572. 

14. Hammer 0, Harper DAT Ryan PD: PAST: Paleontological Statistics 
Software Package for Education and Data Analysis. Palaeontol Electron 
2001,4(1):9-19. 

15. Cadima J, Jolliffe I: On Relationships Between Uncentred And Column- 
Centred Principal Component Analysis. Pak J Statist 2009, 25(4):473-503. 

16. Montenegro Rua EJ, Blanco-Rotea R, Benavides R, Portela C: Santa Eulalia 
De Boveda. Xunta de Galicia;978-84-453-4606-8 2008, 20. 

17. Nomikos P, MacGregor JF: Monitoring batch processes using multiway 
principal component analysis. AlChE J 1994, 40(8):1 361 -1375. 

18. La Gennusa M, Lascari G, Rizzo G, Scaccianoce G: Conflicting needs of the 
thermal indoor environment of museums: In search of a practical 
compromise. J Cult Merit 2008, 9(2):125-134. 



doi:1 0.1 1 86/1 752-1 53X-6-S2-S1 1 

Cite this article as: Visco et al:. Microclimate monitoring in the Career 
Tullianum: temporal and spatial correlation and gradients evidenced by 
multivariate analysis; first campaign. Chemistry Central Journal 2012 6 
(SuppI 2):S11. 



Additional file 1: This contains all the figures with letters as figA 
fig.B cited in the text and necessary to better describe the 
microclimate of the building. 

Additional file 2: Here are stored the raw data obtained from every 
sensor, with also all statistic parameters calculated for every 
column of the main matrix, so for every distribution. In this file 
there are also a copy of the tables cited in the text. The two files 
has been checked, for open, with the major free or open source 
software as MSviewer 2003, LibreOffice, Smartsuite 9.x, with some 
OS. 



